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The properties of the polaron and bipolaron are explored in the ID Jahn- Teller model with dynamical 
quantum phonons. The ground-state properties of the polaron and bipolaron are computed using a 
recently developed variational method. Dynamical properties of the ground state of a polaron are 
investigated by calculating the optical conductivity o{uj). Our numerical results suggest that the 
Jahn- Teller and Holstein polarons are similar. However, in the strong-coupling regime qualitative 
differences in o{u>) between the two models are found and discussed. The influence of the electron- 
phonon coupling and the electrostatic repulsion on the bipolaron binding energy, bipolaron masses, 
and correlation functions is investigated. 
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I. INTRODUCTION 

While research on Jahn- Teller (JT) electron phonon 
coupling spans many decades, a renewed interest in this 
problem has been sparked by the discovery of colossal 
magnetoresistance materials (CMR). In CMR's, it is be- 
lieved, the colossal magnetoresistance effect is a conse- 
quence of the interplay between the double exchange 
mechanism and the lattice effect [1-4]. Furthermore, 
recent studies of alkali doped fullcrcncs A n Ceo indi- 
cate that the Jahn- Teller distortion together with strong 
Coulomb repulsion could be responsible for unusual elec- 
tronic properties of these materials [5-9] . 

Recent advances in computing capabilities have stim- 
ulated development of various numerical techniques such 
as: exact diagonalization techniques (ED) on small 
clusters [10-12], ED on infinite clusters [13,14], quan- 
tum Monte Carlo calculations (QMC) [15,16], variational 
methods [17], and density matrix renormalization group 
techniques (DMRG) [18], that have provided valuable re- 
sults for the Holstein model (HM) in one and more spatial 
dimensions. The most efficient methods provide energies 
for the Holstein polaron problem that are accurate up to 
21 digits in the thermodynamic limit [14]. In addition 
to static quantities, dynamic properties such as spec- 
tral functions and optical conductivity of the Holstein 
model have recently been studied on small lattice clusters 
[19-23]. However, except in the small polaron regime, 
such calculations are subject to pronounced finite-size 
effects. In contrast to the HM, there has been much 
less numerical research devoted to the Jahn- Teller model 
(JTM), mainly due to the far larger Hubert space that 
presents the main obstacle to exact-diagonalization ap- 
proaches. Most of the recent numerical calculations of the 
JTM that take into account the full quantum mechani- 
cal nature of the problem consider only one atom [6,7]. 
A recent path-integral Quantum Monte Carlo approach, 
developed by Kornilovitch, has proven powerful in com- 
puting ground state and spectral properties like the den- 



sity of states of the Jahn- Teller problem [26] . Calculation 
of dynamic properties by that method would require an- 
alytic continuation from imaginary time, which is very 
sensitive to statistical noise. 

The scope of the present work is to compare the static 
properties of the JTM and Holstein model (HM) for the 
case of one and two electrons. In addition, we investigate 
dynamic properties of the Holstein and Jahn- Teller po- 
laron by calculating the optical conductivity. In this case 
we devote equal attention to the Holstein and Jahn- Teller 
models since dynamical properties of the Holstein model 
calculated using an infinite-lattice variational space have 
not been published elsewhere. We use a recently devel- 
oped numerical technique [13,27] to study the JTM for 
the case of one and two electrons on the ID infinite lat- 
tice. Our main goal is to find the numerically exact solu- 
tion of the JTM in the thermodynamic limit. The varia- 
tional method that we use [13,27] is defined on an infinite 
lattice and is not subject to finite-size effects. A standard 
Lanczos method is used to find ground and excited states 
with the selected Hilbert space. The method allows cal- 
culation of physical properties at any wavevector k. In 
the intermediate coupling regime where it is most accu- 
rate, it provides results that are variational in the ther- 
modynamic limit and gives energies accurate to 21 digits 
for the Holstein polaron, up to 7 digits for the Holstein 
bound bipolaron, and 6 digits for the Jahn- Teller polaron. 
While there are no boundary finite-size effects, there are 
nevertheless finite- variational-space effects due to a con- 
straint that only a finite separation between the electron 
and the surrounding phonons is allowed for the polaron, 
and between two electrons and the polaron cloud for the 
bipolaron. Nevertheless, results may be fairly accurate 
even in the two-electron case when electrons are bound 
into a bipolaron, since then electron-electron correlation 
functions decrease exponentially for large separations. 
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II. MODEL 

The simplest, so called E (£> e Jahn- Teller model, con- 
sists of two degenerate electron orbitals and two degen- 
erate phonon modes, 

H = - t J2 (4+i,o, s c i,o,8 + H.c.) + w ]T (atoj + 6+6,) + 
+ gqjp >^[(n^i ?s - n, ;2 ,s)K + + 
(<i, s Ci,2, s + c + i2jS c i>liS )(6 + + 6,)] + 

2,0 

^2 (ni,i,T«i,24 + ni,2,T«i,i,i) > C 1 ) 

where c + Q s creates an electron on site i, orbital o (o = 1,2 
) and spin s, and a + (6 + ) creates dispersionless phonons 
of type a (6) on site i. While the phonon of type a cou- 
ples to the electron density, the mode 6 mediates hop- 
ping between electronic orbitals that are orthogonal in 
the absence of phonon 6. The parameters of the model 
are the intersite hopping matrix element t, dimensionless 
electron-phonon coupling strength g, and optical phonon 
frequency a>o . In the bipolaron case we consider electrons 
with opposite spin. The last two terms in Eq. (1) rep- 
resent the on-site same-orbital (Ui) and different-orbital 
(U2) Coulomb repulsion. For simplicity we have chosen 
U 1 = U 2 = U. 

Most of the calculations will be compared to a simpler 
Holstein model 

H = -t^2 (4+1,8*,* + H - c +UJ oY^ a t a * 

i.s i 

+ SWo^n 1>s (fl+ + <n) + U y^n it fn id . (2) 

i,s i 

While the basic principles of the method have been 
explained elsewhere [13,27] we here present only a brief 
explanation of how the variational space is constructed 
for the case of two electrons in the JTM. Basis states 
for the many-body Hilbert space can be written \<j>) — 
b'1,01, j'2,02; • ■ .nj,n i+ i . ..m„m 1+ i . . .), where the ji, j 2 
and Oi , o 2 are first and second electron site and orbital 
indices respectively, and there are rn phonons of type a 
and rrii phonons of type 6 on site i. A variational sub- 
space is constructed beginning with an initial state where 
both electrons are on the same site and orbital with no 
phonons, and operating repeatedly (Nh times) with the 
off-diagonal pieces (t and g) of the Hamiltonian, Eq. (1). 
The wavefunction is then written in a translation invari- 
ant form. To achieve high accuracy in the strong coupling 
limit where there are many excited phonon quanta, we 
have constructed towers of N a d (N a d — 30 for JTM and 
N a d = 200 for HM) additional phonon excitations (for 



the JT case: one tower for each phonon type) that were 
located on the electron site and on the first neighbor site 
to the left or to the right of the electron position. Such 
towers play a crucial role in achieving convergence in the 
small polaron regime. 

III. JAHN-TELLER POLARON 

A. Static correlation functions 

We start by presenting results of the JT polaron (JTP). 
To investigate the shape of the JTP, we have computed 
various static correlation functions. We first present the 
correlation function for the distribution of the number of 
excited phonons in the vicinity of the electron 

l{i - j) = {i>k\ n t {a+a 3 + 6+6 3 ) \ip k ) , (3) 

where rii — n^i + nj,2 and \ipk) is the polaron wave 
function at momentum k. Numerical results presented 
in Fig. (1) show -f(i — j) calculated at various coupling 
strengths and different wavenumbers. When comparing 
results at different k it becomes evident that the size 
of the polaron at k — tt is larger than that at k = 0. 
As the coupling strength increases, 7(0) increases while 
7(2 ~ J 7^ 0) does not change substantially at k — while 
for larger values of k, -f(i — j 7^ 0) slightly decreases. We 
should stress that due to the the symmetry of the pseu- 
dospin rotation [30], (^ fe |riiO+aj \ip k ) = {^k\n i b + J b.j \ip k ). 
In Fig. (1) we also compare JTP with the Holstein po- 
laron (HP), where 7 is calculated at g — 1.2. While 
differences at small k are very subtle, they become more 
pronounced at larger k where 7(0) of the HP is larger 
than that of the JTP. 
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FIG. 1. j(i — j) vs. (i — j) for the JTP calculated at three 
different coupling strengths and uio = t = 1 (open symbols). 
Results of the HP are presented as black diamonds. 
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Next we present the phonon-phonon correlation func- 
tion of the JTP 



= (i>k\ataibftj \ip k ) 



(4) 



shown in Fig. 2. As expected, phonons of two different 
types are only weakly correlated. At small values of k the 
electron mediated phonon-phonon interaction is always 
attractive (it has a maximum at (i — j) = 0), however, 
at larger values of k and in the weak to intermediate 
coupling regime it becomes repulsive (e(i — j) peaks at 
(i — j) = ±1). In contrast, the expectation for the same 
phonon type, {ip k \ a^aiajaj \tpk), always peaks at (i — 
j) = (not shown). 
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FIG. 2. e(i — j) vs. [i — j) calculated at three different 
coupling strengths and u>o = t = 1. Note that vertical scales 
are different for k = 0, it/4 and k = 7r/2, it. 

The static correlation function between the electron 
position and the oscillator displacement is defined as 



X(i - j) = (ipk\rii(a-j + aj) \tp k ) 



(5) 



We should stress that when using properly sym- 
metrized ground-state wavefunction in the pseudo-spin 
space, we numerically obtain the following equality: 
(ipk\rii(b+ + bj) |V>fc) = (ipk\rii(a+ + aj) \ip k ). In Fig. (3) 
we present \{i — j) in the intermediate coupling regime, 
i.e. g = 1. At k — 0, where the group velocity is zero, 
the deformation around the electron position is limited 
to only a few lattice sites, x(« — j) is always positive 
and seems to be exponentially decaying. At finite but 
small k = it/ 4, the local deformation around the electron 
increases in amplitude. We also notice a ringing effect 
(oscillating deformation) as a consequence of a finite ve- 
locity at k ^ 0. At k = n/2 the ringing effect is strongly 
enhanced while the spatial extent of the deformation in- 
creases in comparison to k = 0. The range of deformation 



reaches its maximum at k — 7r while at the same time 
the lattice displacement on the electron site decreases in 
comparison to smaller values of k. In comparison with 
HP (also presented in Fig. (3)), the main difference is in 
7(0) which in the case of the JTP diminishes strongly 
with increasing k. 
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FIG. 3. x(i - j) vs. (i - j) for JTP and HP, t = 1. 

An important advantage of defining the Hilbert space 
on an infinite lattice is in the ability to calculate the 
energy of the system E at arbitrary value of k. This 
makes our method extremely efficient for computing the 
effective mass of the polaron 



m = 



d 2 E(k ) t 

dk 2 



fe=0- 



(6) 



The second derivative is evaluated by finite differences in 
the neighborhood of k = + dk. Comparison between 
the JTP and HP effective mass is presented in Fig. (4). 
In accordance with the Kornilovitch results [26], the ef- 
fective mass of the HP polaron is smaller than the JTP 
in the weak to intermediate coupling regime. With in- 
creasing coupling the JTP becomes lighter than the HP. 
From Fig. (4) it seems as if the effective mass of the JTP 
and HP displays similar behavior in the strong coupling 
limit. 
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FIG. 4. Effective mass of the JTP and HP vs. g for t = 1, 
u>o — 1 and 5.0, log scale. 

The evolution of a polaron state with changing k can 
be inspected by calculating the quasiparticle weight and 
the mean phonon number Njf 1 : 



^ = Kv fe i4io)i 2 , 



(7) 
(8) 



where c[|0) represents a state with an electron and no 
excited phonons and \ipk) is the solution of the model at 
finite k. Note that the difference between Zk=o and the 
inverse effective mass Eq. (6) comes from k— dependence 
of the polaron self-energy S(fc,w) [14]. At finite k, Zk 
measures the electronic character of the polaronic wave- 
function. In Fig. (5) we present Zk and N% h for the HP 
and JTP calculated at two different coupling strengths 
g = 0.7 and g = 2.0. Consistent with results for the ef- 
fective mass in Fig. (4) , Zk of the JTP is slightly larger 
than the HP at k = 0, g = 2. With increasing k we 
observe a smooth crossover from the predominantly elec- 
tronic character (large Zk) towards the phononic charac- 
ter (small Zk) of a polaron. This crossover is sharper in 
the HP case as Zk of the JTP intersects Zk of the HP 
close to k = 0.3tt (g = 0.7). This intersection is also 
reflected in N% even though N% does not differ sub- 
stantially between the two models. In the case of large 
coupling, g = 2.0, Zk for the JTP is larger than Zk of 
HP which is consistent with the behavior of the effective 
mass in the strong coupling regime. Moreover, 
the JTP is substantially smaller than N% h of HP. 
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FIG. 5. N£ (top) and Z k (bottom) vs. k/n for two differ- 
ent coupling strengths. 



B. Optical conductivity 

Applying standard linear-response theory, the real part 
of the conductivity takes the form (in the limit T = 0), 



Ktr(w) = DS(lo) + a re9 (uj > 0), 



(9) 



where D is the Drude weight at u> = and <r reg is the 
finite-frequency response given by 



E n — E 



S(w + E -E n ), (10) 



where all states \ip n ) for n = 1,2,... belong to k = 
space. It is important to stress that in Eq. (10), the 
current operator J has the same form in the Jahn- Teller 
and Holstcin model, 



J = -ietJ2 ( c t+i,o c i,o -H.c). 



(11) 



Therefore we should expect a similar behavior of a re9 (u>) 
in both models. We next introduce the w-integrated spec- 
tral weight 



S reg (tu) = J a reg {w)du) 
o+ 

and arrive at the sum rule 

S tot - ^f(-H kin ) = S reg + ireHm*- 1 , 



(12) 



(13) 



where S reg = S reg (oo). The sum rule may be used to 
test the numerics by calculating the effective mass by first 
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using Eq. (6) and second, Eq. (13). In all cases presented 
the numerical sum rules are satisfied to at least 6 digits. 

While our numerical method, defined on the infinite 
lattice, enables us to keep the numerical precision in the 
thermodynamic limit of static correlation functions in- 
cluding m* , S tot and S reg within a linewidth, there are 
finite-size effects when presenting finite-frequency quan- 
tities. The reason is that the obtained spectrum is dis- 
crete and we are thus forced to use a small e <~ 0.05 to 
smooth the continued-fraction form of Eq. 10. On the 
other hand, the integrated spectral weight S reg (u}) does 
not depend on e and is thus more reliable. 

We first present <7 reg (ui) and its integral S reg (uj) in 
Fig. (6) for the JT and HM in the weak to interme- 
diate coupling regime. The precision of our results in 
the thermodynamic limit can be tested in the weak cou- 
pling limit, where the optical conductivity threshold is 
known to be at w = w . In this case, other methods, 
defined on finite systems, show pronounced finite-size 
effects due to discreteness in the k — space which are 
reflected in the threshold, larger than uo in the weak- 
coupling limit [19-21], and rather well separated peaks 
corresponding to scattering of the initial k = electron 
state into finite-A: states. The main signatures of a reg (uj) 
at g = 1.0 of both models (JT and HM) are: the spectra 
are strongly asymmetric in frequency, which is a charac- 
teristic of large polarons [24] . Both spectra seem to share 
the same threshold frequency u>q [25] . The first two broad 
peaks are due to one and two phonon emission processes. 
Peaks are more pronounced in the JT than in the HM 
case. With increasing coupling g, (T re9 (uj) of both mod- 
els shifts towards higher frequencies and remains similar. 
The threshold frequency does not change significantly. 

We feel obliged to discuss the constraint of the finite 
Hilbcrt space effects on presented spectra in more detail. 
In the inset of Fig. (6) we present a reg {uo) for the JTP 
calculated for three different sizes of the Hilbert space. 
We can clearly see that broad features (two broad peaks) 
in o~ reg {uj) converge rather well, however, the number 
and position of small peaks that compose broad peaks 
change slightly with increasing N^. Excited states that 
contribute to a reg (w) in the total k = sector can be rep- 
resented as a polaron with the wavevector k p and excited 
phonon or phonons with the total wavevector k p h — — k p . 
Small peaks therefore represent scattering of the polaron 
into different k p states with one or multiple phonon emis- 
sions. Since in our variational space phonons are allowed 
only Nh steps away from the electron, there are only 
discrete k p s allowed for such scattering. In the inset 
of Fig. (6) we present also the integrated conductivity 
S reg (uj) for the same sizes of the Hilbert space. The con- 
vergence in this case is excellent. 

For small hopping t, one can readily calculate optical 
conductivity a re9 {u>) using strong-coupling perturbation 
theory, i.e., the Lang-Firsov transformation. For a two- 
site Holstein model [31], we have 



<7 re9 (uj) = (T n S(uJ — nuj ) 



n=i 



(2g 2 ^ n 



= ne 2 t 2 e- 2g V ^—8{u - ncj ) + 0(t 4 ). (14) 

n— 1 

We note that <j re9 (u) is not quite Poisson distributed. It 
is composed of a series of Dirac 5- functions centered at 
the frequencies u> — nuj + 0(t). Their weights, <r„, are 
second-order in t. For j > 1, the largest a n in Eq. 14 
occurs when n = 2uj (g 2 — 1) ~ 2_E p , which is consistent 
with previous numerical studies [22]. Figure 7a shows 
good agreement between perturbation theory and the nu- 
merical result. 

In contrast, applying perturbation theory to the JTP 
is much more difficult because the single-site problem 
does not have an exact solution. Reasonably accurate 
trial wave functions for the low-lying levels have been 
proposed [32,33]. However, obtaining the analytic wave 
functions for the whole spectrum remains an arduous 
task. In the following, we evaluate the second-order per- 
turbation based on the numerically exact solutions of the 
single-site problem. It is known that the eigenfunctions 
(f> n j of a JT site are characterized by the radial quantum 
number n and the angular quantum number j. Each state 
is doubly degenerate and n = 0, j = ±1/2 for ground 
states. Thus, for a two-site JTP we chose the ground 
state to be 



\lM = ^= [c\\<t>n,l/2) |0, 0) + 4|0, 0) |0 n , 1/2 ) 



(15) 



where |0,0) is the phonon state on the site without the 
electron. We have 



a reg (u>) = weH 2 J2 



\((j) n 1 0,0) (m 1 ,m 2 | (ftp) | 
E n - E Q + (mi + m 2 )u> 



n, mi, m2=0 

x S[u - E n + E - (mi + m 2 )w ] + 0(t 4 ), (16) 

where E n is the eigenvalue of (j> n and n + mi + m 2 > 0. 
The quantum number j is left out because it is conserved 
by current operator. All the nonzero matrix elements ex- 
ist only between states with j = 1/2. In addition, the 
radial excitation of the JTP is anharmonic due to its 
Mexican-hat potential surface. Therefore, the S-function 
peaks are not necessarily located near ruvo- We illustrate 
the distinct feature in Fig. 7b in which the anharmonicity 
manifests itself clearly when the excited states come near 
the cone of the Mexican-hat potential. (The height of the 
cone is g 2 ojQ measured from the bottom of the potential 
surface.) There are actually two peaks near loq and three 
peaks near 2a>o because the excited phonon quanta can ci- 
ther be harmonic or anharmonic. An excellent agreement 
between perturbation theory Eq. 16 and Lanczos diago- 
nalization is found in Fig. 7c. Although the total weight 
of the JTP and HP are approximately equal, in the low- 
frequency regime, we note that S r J^ P (uj) » S t ^p(lo) as a 
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result of the fact that the JTP has a larger quasiparticlc 
weight than the HP [30,33]. 
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FIG. 6. a re9 (uj) and S res {u) vs. w for g = 1 (top) and 
g = 1.5 (bottom) for JTP (full lines) and HP (dashed lines). 
Inset (top) shows a re9 (u>) and 5 re9 (w) for the JTP calculated 
at three sizes of the Hilbert space with Nh = 12, 13, and 14. 
Convergence of S re9 (io) is within the linewidth. 
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FIG. 7. Optical conductivity: second-order perturbation 
theory is compared with the numerical result. The height 
of the peaks represents the weight of the 5-functions. The 
numerical calculations are performed on the 2-site Holstein 
or Jahn- Teller models with t = 0.05 and wo = 1. The elec- 
tron-phonon coupling is g — y/E for all panels. 

Turning to the sum rules presented in Fig. (8), we 
notice a monotonic decrease of the total sum rule 5 tot , 
which indicates a suppression of the electronic kinetic en- 
ergy in the strong coupling regime. The drop in S tot is 
accompanied by a decrease in the Drude weight D, which 
is a measure of the coherent transport properties of a 
polaron. Although the decrease is less steep in the JTP 
case, the two models show strong similarities. The depen- 
dence of the regular part S re9 of the optical conductivity, 
which measures the dissipative part of <j(uj), is not mono- 
tonic in the electron phonon coupling g. It first increases 
with increasing g, reaches the peak near the crossover 
from large to small polaron regimes and then decreases 
with further increase of g. Such an enhancement of S re9 
in the intermediate coupling regime has already been ob- 
served in the Holstein case [19]. In the strong coupling 
limit S tot of both models approach the asymptotic value 
gtot ^ t 2 /u> g 2 (see Fig. (8)), which can be obtained by 
integrating Eq. (14) and summing over all phonon exci- 
tations in the limit g ^> luq. It is rather surprising that 
even though the strong coupling expansion is in principle 
not valid for the JT case, S tot for the JT model seems to 
approach the same asymptotic behavior. All quantities 
in Fig. (8) have an error less than a linewidth. 
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FIG. 8. Partial and total sum rules for wq = 0.5 (top) and 
u)q = 1.0 (bottom) as a function of coupling g, t = 1. Thicker 
lines represent JTP results, thinner lines HP. The thin dotted 
line is a strong coupling result. 
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IV. BIPOLARON 

The following section docs not apply to manganites, 
since the strong Hund's rule coupling renders JT bipo- 
larons unbound. It does, however, apply to other tran- 
sition metal oxides with filled t 2g shells, such as Ni or 
Co. As we will see in this section, similarities between 
the JT and Holstcin models persist in the case of two 
electrons. We start by comparing electron-electron cor- 
relation function 



C(i -j) = (ipk\ riitij \ip k ) 



(17) 



where n>i = n, ; i + n, ; 2 and \tpk) is the bipolaron wave 
function at momentum k. In Fig. (9a) and Fig. (9b) we 
present C(i — j) for k = and different values of the 
electron-phonon and Coulomb interactions. For U = 
the JT bipolaron (JTB) and Holstein bipolaron (HB) 
show very small differences. At smaller coupling g = 1, 
Fig. (9a), the effect of increased U on the JTB and HB 
is similar. The main effect of finite U is to enlarge 
the size of the bipolaron, caused by the electrons try- 
ing to avoid double occupancy. At slightly larger cou- 
pling g = 1.3, Fig. (9b), the effect of finite U is slightly 
more pronounced in the case of JTB than HB. We can 
conclude that the JTB in the intermediate to strong cou- 
pling regime has a slightly larger radius at finite U than 
the HB. 




FIG. 9. C(i - j) vs. i 
and N h = U for JTB. 



• j. We have used N h = 18 for HB 



A larger radius of the bipolaron should imply a larger 
mobility of the bipolaron. It is known that in the HB, 
with increasing U the on-site bipolaron transforms into 
an intra-site bipolaron with a reduced effective mass 
[27-29]. We may therefore expect the JTB to have a 
smaller effective mass than the HB due to a larger radius 



observed in Fig. (9b). This indeed turns out to be the 
case as seen in Fig. (10) where we show the bipolaron 
effective mass vs. U in units of the polaron mass for two 
different coupling strengths. We see that JTB effective 
mass is smaller than HB effective mass by up to a factor 
2 for all U and for both coupling strengths. As shown in 
the inset of Fig. (10), the JTB at g — 1.3 has a smaller 
binding energy than the HB, corresponding to its larger 
radius and smaller mass. The binding energy is defined 
as A = Ebi — 2E po , where Eu and E po are bipolaron and 
polaron ground state energies respectively. 
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FIG. 10. Effective mass of JTB and HB vs. U for lo = 1.0. 
Inset: JTB and HB binding energies vs. U. Finite-Hilbert 
space scaling was used for better accuracy. 

We conclude our numerical investigation of the JTB 
with the phase diagram presented in Fig. (11) for the 
transition from unbound polarons to bound bipolarons. 
We have obtained the phase diagram using the condi- 
tion A(U,g) = 0. It is important to stress that a bound 
bipolaron exists in both models at U = and at any 
finite g. In the weak coupling limit, i.e. j C 1, the 
phase boundary is given by U c — 2cjg 2 7 represented by 
the dashed line in Fig. (11). This result can be derived 
in the classical limit and it is identical for both the JT 
and H models. With increasing g, the JTB phase bound- 
ary deviates slightly downwards from the given analyti- 
cal estimate, while for the HB it deviates upwards. Ap- 
proaching the strong coupling regime, U c (g) for the JTB 
seemingly approaches phase boundary of the HB which 
in the strong coupling regime follows U c = Auig 2 [27]. 
Due to the large Hilbert space of the JTB were unable to 
investigate the strong coupling regime of the JTB more 
precisely in order to determine weather the two phase 
boundaries merge with increasing g. 
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FIG. 11. The phase boundary between two separate po- 
larons and a bipolaron calculated for u>o = 1. 

V. CONCLUSIONS 

We have performed a detailed numerical analysis of 
the JTM in one spatial dimension. The main conclusion 
is that the simple E (g> e JTM is, despite extra orbital 
and phonon degrees of freedom, very similar to the Hol- 
stein model (although the former has a doubly degenerate 
ground state and the latter does not). This result agrees 
with recent quantum Monte Carlo calculations [26] and 
analytical comparison of both models [30] . Closer exami- 
nation reveals that the effective mass of the JTM is larger 
than HP effective mass in the weak coupling regime and 
smaller by roughly a factor of g in the strong coupling 
regime, where the effective masses of JTP and HP share 
the same exponent. 

Turning to spectral properties, our numerical method 
gives extremely reliable results for <j re9 (ui), which can be 
seen in weak to intermediate coupling by observing the 
threshold frequency (the gap in the spectra) . This agrees 
with the weak-coupling prediction, i.e. aj tr — u>o- We 
also support this claim by providing a reg (uj) for systems 
with different numbers of variational states. a re9 (oj) at 
g = 1.0 clearly shows two broad peaks corresponding to 
one- and two- phonon emission. These peaks are much 
more pronounced in the JTM than in the HM. In the 
strong-coupling regime, we find excellent agreement of 
a re9 (cu) for the Holstein model with a simple Lang-Firsov 
strong coupling analytical expression. However, the spec- 
trum of the Jahn- Teller polaron is qualitatively different. 
Due to the anharmonicity of its radial phonon excitation, 
the peaks of a re9 (ui) are not necessarily located near the 
multiples of u>q. For both models, the sum rule in this 
regime is given by S tot <~ t 2 /w g 2 . 



Strong similarities between HM and JTM persist also 
in the case of two bound polarons. Similar to the one po- 
laron case, the Jahn- Teller bipolaron has a smaller effec- 
tive mass than Hubbard bipolaron in intermediate cou- 
pling independent of the strength of U. This is mostly a 
consequence of smaller binding energy. 
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